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Abstract 

We analyze the statistical distribution of neutron stars at the stage of a supersonic propeller. 
An important point of our analysis is allowance for the evolution of the angle of inclination of the 
magnetic axis to the spin axis of the neutron star for the boundary of the transition to the supersonic 
propeller stage for two models: the model with hindered particle escape from the stellar surface and 
the model with free particle escape. As a result, we have shown that a consistent allowance for the 
evolution of the inclination angle in the region of extinct radio pulsars for the two models leads to 
an increase in the total number of neutron stars at the supersonic propeller stage. This increase 
stems from he fact that when allowing for the evolution of the inclination angle x for neutron stars 
in the region of extinct radio pulsars and, hence, for the boundary of the transition to the propeller 
stage, this transition is possible at shorter spin periods (P ~ 5 — 10 s) than assumed in the standard 
model. 



1 Introduction 

Previously, we showed [S] that the transition to the propeller stage could occur at much 

shorter periods than assumed previously when allowance is made for the evolution of the 
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angle of inclination of the magnetic axis to the spin axis of the neutron star. Consequently, 
allowance for the evolution of the inclination angle at the stage of extinct radio pulsars would 
affect significantly the statistics of neuron stars at the propeller stage. 

As and showed, the following two so-called substages that a neutron star passes as 
its spin period decreases must be distinguished: supersonic and subsonic propellers. In both 
cases, the neutron star spins down due to the interaction of the stellar magnetosphere with 
the surrounding matter. 

As a neutron star spins down, the electrodynamic processes cease to play a significant role 
in the pulsar's magnetosphere. When the radius at which the pressure of the ambient medium 
is balanced by the magnetodipole radiation pressure (the Schwarzman radius) becomes equal 
to the Alfven radius, the star passes to the stage of a supersonic propeller. As the spin 
period increases further, the barrier produced by the rotation of the neutron star?s strong 
magnetic field ceases to be efficient, and the plasma begins to penetrate into the pulsar?s 
magnetosphere. Thus, the pulsar passes to the stage of a subsonic propeller. 

Subsequently, due to an even larger spindown, the accretion rate increases and can reach 
its maximum at a certain period P - the neutron star passes to the stage of a steady accre- 
tor. The evolution of a neutron star at the "supersonic propeller-subsonic propeller-steady 
accretor" stages was considered in more detail in the paper by |6j, whose results we used 
here. 

In this paper, we consider the stationary distribution of isolated neutron stars at the 
supersonic propeller stage. As previously jS] , all our calculations were performed for two basic 
models: the model with hindered particle escape from the neutron-star surface (JSl; 0) and 
the model with free particle escape (PP, jTT]). We disregard the possibility of magnetic- field 
evolution and assume that no neutron stars are born at the supersonic propeller stage. An 
important point of our study is that the boundary of the transition of a neutron star from 
the ejector stage to the propeller stage is presented with allowance for the evolution of the 
inclination angle. 

As a result, we obtained the distribution of neutron stars at the supersonic propeller stage 
in spin period P. We show that the number of neutron stars at this stage is larger when 
the evolution of the inclination angle is taken into account for the transition boundary than 
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that in the case where the standard model is considered. We also show once again that when 
the evolution of the inclination angle x is taken into account, neutron stars can pass to the 
propeller stage even at fairly short periods: P ~ 5 — 10 s. 



2 Basic equations 

Lipunov and Popov [H] formulated an important assertion: for a constant magnetic field, 
the ejector stage at reasonable parameters is always longer than the propeller stage. Our 
objective in this part of the paper is to determine the number of neutron stars at the propeller 
stage disregarding both the magnetic-field decay and the evolution of the angle of inclination 
of the magnetic axis to the spin axis of the star at this stage. 

A neutron star at the propeller stage is known to spin down due to the transfer of angular 
momentum to the surrounding matter ([^, [7]). There are many formulas that describe this 
spindown (see |H], [Oj)- Nevertheless, virtually all of them can be reduced to the form [12] 

diu u2 

Here / ~ 2/5Mi?^ is the moment of inertia of the star, /x is the magnetic dipole moment, 
fct ~ 1 is a dimensionless constant, and Ra is the Alfven radius: 
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where M is the mass of the neutron star. The acceleration rate can be estimated as 

M ~ TiRlpV^, (3) 

where Voo is the space velocity of the neutron star relative to the surrounding interstellar 
medium, p is the density of the interstellar medium, and Ra is the gravitational capture 
radius, which for spherically symmetric accretion is 

= -r?^- (4) 

QO 

The maximum possible accretion rate onto an isolated neutron star (for spherically sym- 
metric accretion) was determined by Ikhsanov [U] and can be specified by the following 
expression: 

M < M^ax ^ 2.3 ■ 10"p_24mV6-3gs-\ (5) 
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Figure 1: Evolution of neutron stars under the assumption of hindered particle escape: I - active 
radio pulsars, II - extinct radio pulsars, and III - pulsars at the supersonic propeller stage. 
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Figure 2: Evolution of neutron stars under the assumption of free particle escape. The notation is 
the same as that in Fig. ^ 

where m is the mass of the neutron star in terms of solar masses, p_24 = p/10~^^ g cm~^, and 
Vq = \4o/10® cm s~^. Given all of the factors listed above, the transfer of angular momentum 
to the surrounding interstellar medium is found to cause the neutron star to spin down at 
the rate 

^ = 7.3 • lO-'^Bli'P^ (6) 

Further, it should be noted that for the model with hindered particle escape from the 
pulsar's surface, the transition to the supersonic propeller stage for most neutron stars can 
occur only at fairly small angles x between the magnetic axis and the spin axis (Fig. Con- 
sequently, in the case of hindered particle escape, the neutron stars passing to the propeller 
stage are almost aligned rotators, and the subsequent evolution of the inclination angle is 
unimportant. As a result, it can be assumed with a sufficient accuracy for the Ruderman- 
Sutherland model that the inclination angle does not evolve at the supersonic propeller stage. 

As regards the model with free particle escape from the neutron-star surface, the charac- 
teristic period at which the transition to the propeller stage occurs is much longer than the 
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death period. Consequently, most of the neutron stars will pass to the region of angles close 
to 90° before the transition of the star to the supersonic propeller stage (see Fig. |2I). 

Thus, we are interested in the statistics of pulsars at the supersonic propeller stage without 
considering the angle of axial inclination. Accordingly, to determine the stationary distribu- 
tion function of neutron stars, let us introduce their distribution N3{P,B) in period P and 
magnetic field B. Next, as in the case of extinct radio pulsars, we assume that the function 
N^lPjB) depends on the magnetic field B of the neutron star only as a parameter; i.e., we 
disregard the magnetic-field evolution over the lifetime of the neutron star. As a result, the 
kinetic equation for neutron stars at the supersonic propeller stage can be written as 

where the source function F{P, B) is determined by the neutron-star flux from the region of 
extinct radio pulsars. 

Finally, we must specify the boundary periods for the supersonic propeller stage, i.e., the 
period at which the transition from the ejector stage to the propeller stage occurs and the 
period at which the accretion of matter from the interstellar medium onto the neutron star 
begins. As we found previously [3], the transition from the ejector stage to the propeller 
stage occurs for the model with hindered particle escape from the neutron-star surface at the 
spin period 

Ppr = PESin^'\. (8) 
At the same time, the following expression holds for the model with free particle escape: 

Ppr = PeCOS^'\. (9) 

Below, the expression for the limiting period Pe has the standard form ^U] 

^^RfRc^V'^v^ ( Bl y/\^ 

C \GM) C y^TiPooVl,) 

1/2 1/2 1/2 1/2 

^ \ 1 n2 A'-SO ^7 V-j 



x(-) ^10- 7- s. (10) 



Here, poo and c^o are, respectively, the density and the speed of sound in the interstellar 
medium (cy and v^ are in the units of 10^ cm s^^); B1^^ = SvrpooC^ is the corresponding 
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magnetic energy density ((-Bcxt)-6 is in units of 10~^ G); and /lao is the magnetic moment of 
the neutron star in units of 10^'' G cm^. In addition, we assumed that the velocity of neutron 
stars foo is much larger than the speed of sound in the interstellar medium. 

Next, it is necessary to determine the spin period at which the star passes to the subsonic 
propeller stage. For this purpose, we recall that a pulsar can be at the supersonic propeller 
stage as long as its angular velocity is large enough, and, accordingly, the centrifugal ac- 
celeration exceeds the free-fall acceleration. Consequently, the condition under which the 
accretion of matter from the interstellar medium onto an isolated neutron star begins can be 
specified as follows: 



u:^Ra = (11) 



GM 

where = 27r/P is the angular velocity of the neutron star. 

As a result, for the boundary period corresponding to the transition of an isolated neutron 
star to the subsonic propeller stage, we obtain an expression that is similar to that derived 
by Ikhsanov [Hj: 

Pe, ^ 2A^i%'M^^"m''/' s, (12) 

where M15 is the accretion rate in units of 10^^ g , and m = M/Mq. 

Substituting the parameters into the above formula, we can estimate the period at which 
the neutron star leaves the supersonic propeller stage: 

Peg ^ erOPi'f S. (13) 

Let us now consider the spin period distribution of neutron stars passed to the supersonic 
propeller stage for the standard model, i.e., for the model in which the evolution of the 
inclination angle in the region of extinct radio pulsars is disregarded. In this case, the kinetic 
equation can be written as 

d ( dP\ , , 

It should also be recalled that N2{P) = NfP for the region of extinct radio pulsars in 
the standard model, while the derivative of the spin period at the supersonic propeller stage 
can be represented by Eq. (jH)). As a result, we obtain the following distribution of neutron 
stars in spin period P for the model where the evolution of the angle of inclination of the 
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magnetic axis to the spin axis is disregarded: 

Ns{P) = 1.4 ■ 10-'^NfP~\ if P > Pe, 
N3{P) = 0, iiP<PE, (15) 

In this case, the total number of neutron stars at the supersonic propeller stage for the 
model in which the inclination angle is disregarded can be estimated a.s = 1 .4 x 10~^ A^j / Pg. 

3 Model with hindered particle escape from the neutron- 
star surface 

Let us consider the statistics of pulsars at the supersonic propeller stage in more detail 
assuming that the model with hindered particle escape from the neutron-star surface is valid. 
In this case, the boundary of the transition from the region of extinct radio pulsars (the 
ejector stage) to the supersonic propeller region is defined by Eq. (jHJ. 

Thus, the source function F{P, B) in the kinetic equation ((Tj) can be written as 



F(P,P) = F(P,i?,x)L,„,=p.p-., (16) 

where 

F{P,B,x)=xN2. (17) 

The distribution function N2{P, B,x) in the latter expression refers to the region of extinct 
radio pulsars, and we found it previously Recall that this function for the model with 
hindered particle escape from the neutron-star surface can be specified as 

N,[P,B,x]=0.7k^NfB%\l + Bu)-'-'^^^^^^, (18) 

sznx 



where 



GiO=^-^i^-5+"^M\0- (19) 



In this case, we used two auxiliary functions: 



c(0 = B^'^e-' = cosx, 
s(0 = (1 - c\OY^^ = smx 
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To determine the derivative of the angle of inchnation of the magnetic axis to the spin 
axis of the neutron star, we use the integral of motion in explicit form. If the Ruderman- 
Sutherland model is valid, then the plasma in neutron stars at the stage of extinct radio 
pulsar will fill only the inner regions of the magnetosphere. Consequently, the energy losses 
for such pulsars can be assumed to be identical to the magnetosphere losses. Thus, the 
quantity / = cosx/ P is conserved during the evolution. As a result, we obtain the derivative 
of the angle of axial inclination in the region of extinct radio pulsars 

X = -^cotx. (20) 

The differential equation for the distribution function N3{P,B) of the pulsars at the 
supersonic propeller stage then takes the form 

A (iV3P-) = OMk^N.BTil + B,,r'-j,G (I) . (21) 

Here, we denoted y = (1 - P^Pe'^Y^^. 

The sought-for distribution function can then be represented as 

p 

Ns{P, B) = OMkNNfBlfil + Bu)-'-'P-' J dP, (22) 

P 

where /cat ~ 4.4 is the normalization factor, and Nj is the number of normal radio pulsars 
far from the death line (for more detail, see [3|). 
Since the death line of normal radio pulsars 

Pd{x) = BT{cosxr' s (23) 

crosses the boundary of the transition to the supersonic propeller stage at a fairly small angle 
X {sinx = P^^P°'^^^), the point of intersection of the death line for normal radio pulsars (j^Hjl 
with the sinx = axis can be taken with a good accuracy as the point with the minimum 
period Pmin (see Fig. [T]). It thus follows that 

Pm^n = B%'' S. (24) 

Since the period at which a pulsar passes from the ejector stage to the propeller stage is 
defined by Eq. (jSJ and the period at which a pulsar leaves the supersonic propeller stage is 
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Figure 3: Distribution of neutron stars at the supersonic propeller stage N^{P) for the model with 
hindered particle escape. 

defined by Eq. ()13|) . we derive the following expression for the period distribution function 
of neutron stars at the supersonic propeller stage: 



pl5/8 



B 



N:i(P,B) = Omk^NfP J 

5-10-4p7/6 



20/7, 



1 + 5i2)-'-^ X 



(25) 
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Figure 01 shows the distribution of pulsars at the supersonic propeller stage in period P 
for three different values of Pe- 10, 30, and 100 s. This distribution was constructed by 
assuming that the model with hindered particle escape from the neutron-star surface was 
valid. In addition, we assumed that no neutron stars are born at the supersonic propeller 
stage and did not consider the possibility of magnetic-field evolution over the lifetime of the 
pulsar at this stage. As we see from this figure, the number of such neutron stars is much 
smaller than the number of both active and extinct radio pulsars. 

As was noted above, when allowance is made for the evolution of the angle of inclination 
of the magnetic axis to the spin axis in the region of extinct radio pulsars (and, accordingly, 
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for the boundary of the transition to the supersonic propeller stage), the transition to the 
propeller stage can occur even at fairly short spin periods of the neutron star. The numerically 
computed distribution of pulsars at the supersonic propeller stage shown in the figure confirms 
this conclusion. We also see from Fig. El that the total number of neutron stars increases 
compared to the standard model when allowance is made for the evolution of the inclination 
angle in the region of extinct radio pulsars and, accordingly, for the boundary of the transition 
to the propeller stage. 

It should also be noted that if the model with hindered particle escape is valid, the 
limiting period Pe affects only slightly the distribution of neutron stars at the supersonic 
propeller stage. This situation arises from the fact that most of the neutron stars crossing 
the boundary of the region of extinct radio pulsars sinx = / have short spin periods, 
P < 100 s (Fig.H}. 

4 Model with free particle escape from the neutron- 
star surface 

Recall that the boundary of the transition from the region of extinct radio pulsars to the 
supersonic propeller stage for the model with free particle escape from the neutron-star 
surface is specified by Eq. (see Fig. |2I)- Consequently, the following relation is valid for 
the source function F[P, B) determined by the neutron-star flux from the region of extinct 
radio pulsars: 

F{P,B)=F{P,B,x)los^=p.p-^. (26) 

Therefore, as for the Ruderman-Sutherland model, the source function F[P, B) can be written 
as 

F{P,B,x) = xN2. (27) 

If the Arons's model is valid for the existing neutron stars, then, as was mentioned above, 
the neutron-star magnetosphere for extinct radio pulsars still remains completely filled with 
plasma, which screens the magnetic field. As a result, all of the energy losses are attributable 
to longitudinal currents, while the magnetodipole losses are completely screened. In this case. 
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the quantity Id = sin xl P is conserved during the evolution at the ejector stage (including the 
region of extinct radio pulsars). It thus follows that the derivative of the angle of inclination 
of the magnetic axis to the spin axis of the neutron star is 

P , , 

X = -ptanx- (28) 

Based on our previous results, we use the following expression for the distribution function 
of extinct radio pulsars 

N,{P, B, X) = O.Sk^Nj , ^^'^'^^+^^^^ (29) 

cosx (l + O.lBli {sinx/P) ) 

Let us change the variable 

Substituting (|^. ^I7\ . (j2HI), and into the kinetic equation ((7j) yields the differential 
equation for N^i^P, B) 

^(N^P^) = 4 X IQ-^ PskNNfBf ^.7 A{1 + Bi2Y^-' p-^"^^ x 

y 



{i+o.iBW{y/pry"' 

The solution of this equation is the sought-for distribution function 

NsiP, B)=4x 10-^ PEkNNfBl^'\l + B12Y''' P'^ x 
p 

f y dP 
\L ^(i + o.iB!,/fe/P)")"-'- 

To find minimum spin period of a neutron star that passes to the supersonic propeller 
stage, let us turn to Fig. El We see from this figure that the point of intersection between 
the line of transition of pulsars to the region where the angle of inclination of the magnetic 

1 /2 

axis to the spin axis is close to 90° {cosx = {^R/c) ) and the line of transition from the 
region of extinct radio pulsars to the propeller stage Q, where Vt = 27r/P is an angular 
velocity, corresponds to this period. Thus, the minimum spin period of a neutron star at the 
supersonic propeller stage is 

P^,„ = 0.18AA'/''B%\ (32) 
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Figure 4: Distribution of neutron stars at the supersonic propeller stage N3{P) for the model with 
free particle escape. 

1/2 

Here, we represented the limiting period Pe as Pe = ■ 

As a result, we obtain the spin period distribution of neutron stars at the supersonic 
propeller stage for the model with free particle escape from the pulsar's surface by taking 
into account the boundaries of the transition of the neutron star to the supersonic propeller 
stage (jni) and the subsonic propeller stage (fT3|l and using the distribution function (|3T|l : 

69A-2p5/2 

-3d U at D-2 f D2.74/1 , D ^-3.7 



N,{P) = 4 X W-'PEk^NfP-^ J Bl^\l + B,2) - X 

y dP 



0.185P0» 



P2^/i^(l + 0.15iY(l//P)3-i 



>0.9 • 



Recall that in this expression, the normalization factor is ki^ 4.4, and Nf is the number of 
normal radio pulsars far from the death line. 

Figure |3] shows the distribution of neutron stars at the supersonic propeller stage in spin 
period P under the assumption that the Arons's model is valid for three limiting periods 
Pe'- 10, 30, and 100 s. In addition. Figs. El and El separately show the distributions for 
Pe = 30 s and Pe = 100 s, respectively. As one can see from these figures, when allowance 
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Figure 5: Distribution of neutron stars at the supersonic propeller stage N^{P) for the model with 
free particle escape for the limiting period Pe = 30 s. 
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Figure 6: Same as Fig. El for Pe = 100 s. 
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is made for the evolution of the inchnation angle x the ejector stage and, accordingly, 
for the boundary of the transition to the propeller stage, the number of neutron stars at 
the supersonic propeller stage is indeed found to be significant at fairly short spin periods 
(P ~ 5 — 10 s). Thus, a consistent allowance for the angle of inclination of the magnetic axis 
to the spin axis of the neutron star at the ejector stage leads to an increase in the number of 
pulsars that passed to the supersonic propeller stage. 

It should also noted that, in contrast to the model with hindered particle escape from 
the neutron-star surface, the limiting period Pe plays a significant role in the distribution 
of neutron stars at the supersonic propeller stage for the Arons's model. This is because 
the distribution of neutron stars at the boundary of the transition from the region of extinct 
radio pulsars to the propeller stage is virtually uniform for this model. 

5 Conclusions 

We have determined the distributions of neutron stars at the supersonic propeller stage in 
spin period P taking into account the fact that a consistent allowance for the evolution of 
the angle of inclination of the magnetic axis to the spn axis for the region of extinct radio 
pulsars directly affects the statistics of neutron stars at the supersonic propeller stage. All 
our calculations were performed for two models of the particle acceleration region: the model 
with hindered particle escape |T^ and the model with free particle escape p. 

As a result we have made sure once again that a consistent allowance for the evolution of 
the inclination angle x the region of extinct radio pulsars both for the model with hindered 
particle escape and for the model with free particle escape leads to an increase in the total 
number of neutron stars at the supersonic propeller stage. This increase stems from the fact 
that when the inclination angle x fo^^ the evolution of neutron stars in the region of extinct 
radio pulsars and, hence, for the boundary of the transition to the supersonic propeller stage 
is included, this transition is possible at shorter spin periods of the neutron star (P ~ 5 — 10 s) 
than assumed in the standard model. 
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